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own  manuals  or  on-line  documentation  system.  In  order  to  determine  what  software  is  avail- 
able to  solve  a particular  problem,  users  must  search  through  a very  large,  heterogeneous 
collection  of  information.  This  is  a tedious  and  error-prone  process.  As  a result,  there  has 
been  much  interest  in  the  development  of  automated  advisory  systems  to  help  users  select 
software. 

Keyword  search  is  a popular  technique  used  for  this  purpose.  In  such  a system  keywords 
or  phrases  are  assigned  to  each  piece  of  software  to  succinctly  define  its  purpose,  and  the 
set  of  aU  such  keywords  axe  entered  into  a database.  Keyword-based  selection  systems 
query  users  for  a set  of  keywords  and  then  present  a fist  of  software  modules  which  contain 
them.  A major  difficulty  with  such  systems  is  that  users  often  have  trouble  in  providing  the 
appropriate  keywords  for  a given  mathematical  or  statistical  problem.  There  is  such  a wealth 
of  alternate  mathematical  and  statistical  terminology  that  it  would  be  a rare  occurrence  for 
two  separate  knowledgeable  persons  to  assign  the  same  set  of  keywords  to  a given  software 
module.  Users  of  these  systems,  who  are  usually  much  less  familiar  with  the  terminology, 
often  find  that  there  are  either  too  many  software  modules  associated  with  the  keywords  that 
they  have  specified,  or  none.  One  can  attempt  to  ameliorate  this  problem  by  imposing  a 
standard  set  of  keywords  or  by  implementing  a very  elaborate  keyword-specification  scheme 
for  users.  The  former  is  difficult  to  maintain  and  the  latter  may  not  be  easy  to  use. 

Classification  systems  have  long  been  used  to  give  structure  to  large  bodies  of  infor- 
mation. A well-formulated  system  can  improve  understanding  of  the  information  as  well  as 
ease  access  to  it,  thus  making  the  information  more  useful.  The  Dewey  Decimal  System, 
for  example,  provides  a means  for  librarians  to  maintain  a large  collection  of  books.  Since 
the  system  is  subject-oriented,  hbrary  users  can  quickly  find  books  in  a given  subject  area. 
Likewise,  a subject-oriented  classification  system  can  be  an  effective  means  of  directing  users 
to  appropriate  mathematical  and  statistical  software. 

To  be  effective,  such  a classification  system  must  have  the  following  properties. 

1.  Problem-orientation.  It  must  classify  the  problems  which  can  be  solved  by  computer 
software.  Other  orientations,  such  as  classification  by  algorithm  or  classification  by 
software  package,  are  of  less  interest  to  end  users. 

2.  Variable- level  tree  structure.  A tree  structure  is  the  most  natural  for  a classification 
system.  Allowing  arbitrary  levels  of  refinement  permits  the  system  to  adapt  to  both 
mature  and  young  subject  areas.  In  young  subject  areas  little  software  is  available,  and 
hence  little  refinement  is  necessary.  In  mature  areas  where  much  software  is  available, 
increased  refinement  is  necessary  to  distinguish  among  the  choices. 

3.  Active  maintenance.  The  system  must  be  monitored  and  revised  over  time  to  best 
reflect  the  current  state  of  the  rapidly  growing  mathematical  and  statistical  software 
collection.  The  tree  structure  facihtates  this  by  insuring  that  modifications  of  the 
system  are  localized. 

The  classification  system  can  be  thought  of  as  a pre- defined  hierarchy  of  keywords.  Since 
the  entire  universe  of  these  keywords  is  visible  to  both  developers  and  users  the  possibility 
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of  finding  the  correct  software  is  improved.  In  addition,  the  tree  structure  shows  explicit 
relationships  among  the  set  of  keywords  which  can  aid  in  users’  understanding  of  them. 

To  use  such  a system,  each  piece  of  software  must  be  assigned  a class.  Classification- 
based  software  selection  systems  permit  users  to  incrementally  refine  their  specification  of 
the  problem  by  using  the  classification  system  as  a decision  tree.  When  the  most  appropriate 
problem  class  has  been  selected  by  the  user,  software  which  contains  the  selected  classification 
is  presented. 

Classification  systems  are  not  themselves  free  of  problems.  For  example,  they  partition 
software  much  more  coarsely  than  do  keywords,  and  users  may  not  be  fluent  in  the  termi- 
nology of  the  classification  system.  Nevertheless,  we  beheve  that  the  structure  imposed  by 
such  systems  improves  both  the  user’s  access  to  information  and  the  implementor’s  job  of 
maintaining  it. 

In  this  paper  we  describe  a particular  classification  system  for  mathematical  and  sta- 
tistical software  which  meets  the  criteria  described  above.  The  system  is  an  outgrowth  of 
the  Guide  to  Available  Mathematical  Software  (GAMS)  project  at  NIST.  We  begin  by  de- 
scribing the  origins  of  system,  and  then  outline  its  current  version.  We  next  discuss  how 
such  a system  can  be  effectively  used  for  software  selection.  The  entire  classification  system 
is  included  in  Appendix  A.  Appendix  B is  devoted  to  summarizing  the  differences  between 
the  current  version  and  its  last  widely-publicized  version. 


2 Origins  of  the  System 

The  system  described  in  this  paper  has  its  origins  in  a software  classification  scheme  devised 
in  the  1960s  by  SHARE,  the  IBM  Users  Group.  The  SHARE  scheme  had  a fixed  two- 
level  structure  which  led  to  very  broad  classes.  Such  a coarse  structure  was  not  suited  to 
classifying  very  large  software  collections. 

In  1975  John  Bolstad  proposed  a substantial  revision  of  the  SHARE  scheme  which 
eliminated  many  of  its  weaknesses  [7].  The  Bolstad  system  was  a multi-level  tree-structured 
scheme.  Unfortunately,  it  attempted  to  maintain  compatibihty  with  the  SHARE  scheme 
wherever  possible  and  hence  it  inherited  much  of  the  SHARE  system’s  illogical  organization. 
Nevertheless,  the  scheme  was  a great  improvement  and  versions  of  it  were  adopted  for  use 
by  a number  of  groups,  including  NBS^. 

Many  difficulties  in  using  the  Bolstad  scheme  surfaced  at  NBS  while  attempting  to 
classify  about  2500  subprograms,  representing  most  of  the  widely-distributed  mathematical 
software  then  available.  As  a result,  many  new  classifications  were  added  and  some  sections 
were  completely  reworked.  The  resulting  scheme  was  used  in  the  first  GAMS  software  catalog 
[5]  and  was  adopted  for  use  in  the  documentation  of  the  SLATEC^  Common  Math  Library 
[10]. 

In  1983  Boisvert,  Howe  and  Kahaner  completely  revised  the  classification  system  used 
in  GAMS  and  published  it  for  public  review  [2].  This  became  known  as  version  1.0  of  the 
GAMS  Classification  Scheme  for  mathematical  and  statistical  software.  The  new  system  kept 

^National  Bureau  of  Standards.  NBS  became  NIST  in  1988. 

^Sandia-Los  Alamos-Air  Force  Weapons  Laboratory  Technical  Exchange  Committee 
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the  Bolstad  philosophy  while  providing  an  organization  which  more  accurately  reflected  the 
then  current  state  of  mathematical  and  statistical  software.  In  addition,  many  sections  of 
the  Bolstad  scheme  not  directly  related  to  mathematical  or  statistical  software  were  deleted. 
Instead,  the  new  scheme  was  viewed  as  a node  in  a larger  scheme  which  encompassed  all 
computer  software. 

Minor  modiflcations  of  the  GAMS  scheme  have  appeared  in  [3]  (version  1.1),  [4]  (version 
1.2),  and  [8]  (version  1.3).  Since  1983,  the  scheme  has  been  adopted  for  use  by  a number 
of  institutions,  including  Amoco  Production  Research,  ASA  Statistical  Computing  Section 
Committee  on  Statistical  Algorithms,  C.  Abaci  (The  Scientific  Desk),  Centro  di  Calccolo 
Elettronico  Interuniversitario  deUTtalia  Nord-orientale  (CINECA,  Bologna),  Eigen bssische 
Technische  Hochschule  Zurich  (ETH)  Seminar  fur  Angewandte  Mathematik,  IMSL  Inc., 
Konrad- Zuse-Zentrum  fiir  Informationstechnik  Berlin  (ZIB),  Los  Alamos  National  Labora- 
tories, National  Center  for  Atmospheric  Research,  SLATEC  Common  Math  Library  Subcom- 
mittee, Stanford  Linear  Accelerator  Center  (SLAC),  State  University  of  Utrecht  Academic 
Computer  Center,  and  University  of  Texas  System  Center  for  High  Performance  Computing. 

3 Version  2 of  the  Classification  System 

Much  new  mathematical  and  statistical  software  has  appeared  since  the  GAMS  classifica- 
tion system  was  developed.  In  many  cases  this  software  addresses  problems  not  exphcitly 
included  in  the  original  classification  system.  In  other  cases,  significant  new  software  pack- 
ages have  appeared  which  have  provided  improved  methods  of  organizing  certain  subject 
areas.  Because  of  this  we  have  again  found  it  necessary  to  modify  the  GAMS  Classifica- 
tion System.  The  resulting  scheme,  termed  version  2.0,  appeared  in  the  recently  published 
GAMS  software  catalog  [6],  and  is  reproduced  in  its  entirety  in  Appendix  A of  this  paper. 
Appendix  B outfines  the  differences  between  this  version  and  Version  1.2. 

The  highest  levels  of  the  classification  system  have  remain  unchanged  since  version  1.0. 
In  the  following  we  describe  the  purpose  of  each. 

• A.  Arithmetic,  error  analysis 

Contains  software  implementing  elementary  arithmetic  operations  on  non-standard 
data  types.  Examples  are  extended  precision  arithmetic  and  interval  arithmetic.  Also 
included  are  systems  or  utilities  which  do  general-purpose  error  analysis.  Finally, 
software  for  accelerating  the  convergence  of  sequences  is  also  found  here. 

• B.  Number  theory 

Software  classified  here  performs  such  number-theoretic  calculations  as  the  decompo- 
sition of  integers  into  prime  factors. 

• C.  Elementary  and  special  functions 

Software  for  evaluating  both  elementary  and  specialized  mathematical  functions  is 
found  here.  Examples  of  elementary  functions  are  trigonometric  functions,  exponen- 
tials, and  polynomials.  Examples  of  special  functions  are  Bessel  functions  and  Gamma 
functions.  Statistical  functions  such  as  probability  density  functions  are  found  in  class 
L5. 
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• D.  Linear  Algebra 

This  class  includes  elementary  vector  and  matrix  operations,  matrix  factorizations, 
solution  of  hnear  systems,  eigenvalue  problems,  determinants,  and  inverses. 

• E.  Interpolation 

Software  for  finding  a function  which  “passes  through”  given  data  values  in  one  or 
more  dimensions  is  found  here.  K the  data  have  noise  then  classes  K or  L8  are  more 
appropriate. 

• F.  Solution  of  nonlinear  equations 

This  class  contains  software  for  solving  systems  of  nonlinear  equations.  Software  for 
single  nonlinear  equations  and  polynomial  equations  are  also  included. 

• G.  Optimization 

Software  for  minimizing  or  maximizing  functions  with  or  without  constraints  is  found 
here.  This  includes  linear  programming,  nonhnear  programming,  integer  programming, 
network  optimization,  and  optimal  control. 

• H.  Differentiation,  integration 

Here  one  finds  software  for  estimating  derivatives  and  evaluating  integrals. 

• I.  Differential  and  integral  equations 

Software  for  solving  ordinary  differential  equations,  partial  differential  equations,  and 
integral  equations  is  found  here. 

• J.  Integral  transforms 

This  class  includes  software  for  Fourier  transforms,  trigonometric  transforms,  Laplace 
transforms,  Hilbert  transforms,  convolutions,  etc. 

• K.  Approximation 

Software  for  determining  best  approximations  to  functions  or  data  in  various  norms 
(e.g.,  Z/i,  L2,  Loo)  are  classified  here.  Software  for  approximation  followed  by  statistical 
analysis  (e.g.,  regression)  is  classified  in  L8.  Software  for  solving  linear  algebraic  sys- 
tems (e.g.  solution  of  overdetermined  systems  in  the  least  squares  sense)  is  classified 
in  D9.  Software  for  interpolation  is  classified  in  E. 

• L.  Statistics,  probability 

Software  for  statistical  computing  is  classified  here.  This  includes  data  summarization 
and  manipulation,  elementary  data  analysis  (e.g.,  calculating  the  sample  mean),  sta- 
tistical graphics,  statistical  function  evaluation,  random  number  generation,  analysis 
of  variance,  regression,  categorical  data  analysis,  time  series  analysis,  correlation  anal- 
ysis, discriminant  analysis,  covariance  structure  models,  cluster  analysis,  and  survival 
analysis. 

• M.  Simulation,  stochastic  modeling 

Software  for  building  and  studying  stochastic  models  is  classified  here. 
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• N.  Data  handling 

Data  handling  includes  various  operations  such  as  input,  output,  sorting,  searching, 
merging,  and  permuting.  Software  implementing  useful  data  structures  such  as  heaps 
and  trees  is  also  found  here. 

• 0.  Symbolic  computation 

Software  for  manipulating  mathematical  expressions  in  their  symbolic  form  is  classified 
here. 

• P.  Computational  geometry 

This  class  includes  software  for  fundamental  geometric  calculations  (e.g.,  areas  and  vol- 
umes) and  implementation  of  algorithms  for  geometric  problems,  such  as  computation 
of  the  convex  hull  and  the  Voronoi  diagram. 

• Q.  Graphics 

General-purpose  computer  graphics  is  classified  here.  Statistical  graphics  is  in  L3. 

• R.  Service  routines 

This  class  includes  software  which  performs  low-level  utihty  functions  such  as  error 
checking,  error  handling,  and  retrieval  of  information  about  machine  characteristics. 

• 5.  Software  development  tools 

Tools  which  facihtate  mathematical  software  development  and  maintenance  and  clas- 
sified here.  Tool  types  include  program  transformation  (e.g.,  convert  to  double  pre- 
cision), static  analysis  (e.g.,  flow  analysis,  interface  analysis),  and  dynamic  analysis 
(e.g.,  tracing,  timing,  assertion  checking). 

• Z.  Other 

This  class  contains  software  which  does  not  fit  anywhere  else. 

Strictly  speaking,  classes  N,  Q,  R,  and  S are  not  mathematics  or  statistics.  They  have 
been  included  because  such  software  is  commonly  found  in  mathematical  and  statistical 
packages. 


4 Using  the  System 

Version  2.0  of  the  GAMS  Classification  System  has  been  used  to  classify  more  than  5000 
software  modules  from  40  separate  packages  at  NIST.  The  1990  edition  of  the  Guide  to 
Available  Mathematical  Software  [6]  lists  all  of  these  modules  in  order  of  their  classifications. 
A software  advisory  system  called  the  GAMS  Interactive  Consultant  (GAMSIC)  provides  on- 
line access  to  this  same  information  [1].  Such  successful  application  indicates  that  the  system 
adequately  reflects  the  current  state  of  mathematical  and  statistical  software  development. 
Some  of  the  well-known  packages  which  have  been  classified  are  fisted  in  Table  1.  In  what 
follows  we  discuss  some  of  the  guidelines  which  have  governed  the  use  of  the  system  at  NIST. 
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Table  1:  Some  well-known  packages  classified  at  NIST 


Package 

n^ 

Description 

BMDP 

43 

Programs  for  statistical  data  analysis. 

CALGO 

172 

The  Collected  ALGOrithms  of  the  ACM.  Programs  published  by 
the  ACM  Transactions  on  Mathematical  Software.  (1975-88) 

CMLIB 

763 

The  NIST  Core  Math  Library.  A collection  of  public-domain 
Fortran  subroutine  packages.  (Includes  LINPACK,  EISPACK, 
FISHPAK,  QUADPACK,  FFTPKG,etc.,  many  of  which  are  also 
found  in  the  SLATEC  hbrary) 

DATAPAC 

169 

Fortran  subprograms  for  statistical  data  analysis. 

Dataplot 

87 

Interactive  graphical  and  statistical  data  analysis  program. 

IMSL 

470 

Fortran  subprograms  for  mathematics  and  statistics.  (Version 
9.2) 

MATH/LIBRARY 

668 

Fortran  subprograms  for  mathematics  from  IMSL  Inc.  (Version 
1.0) 

STAT/ LIBRARY 

620 

Fortran  subprograms  for  statistics  from  IMSL  Inc.  (Version  1.0) 

SFUN/LIBRARY 

297 

Fortran  subprograms  for  evaluating  special  functions  from  IMSL 
Inc.  (Version  2.0) 

MAGEV 

80 

The  MAth/GEophysical  Vector  library.  Fortran  subprograms  de- 
veloped for  the  Cyber  205.  (Version  3.3) 

MINITAB 

63 

Program  for  statistical  data  analysis. 

NAG 

774 

Fortran  subprograms  for  mathematics  and  statistics.  (Mark  13) 

NMS 

50 

Fortran  subprograms  from  [9]. 

PORT 

327 

Fortran  subprograms  for  mathematics.  (Version  2) 

SAS 

40 

Program  for  statistical  data  analyses.  (Version  15.8) 

Scientific  Desk 

329 

Fortran  subprograms  for  mathematics  and  statistics  for  use  on 
PC’s.  (Version  4) 

SPSS 

28 

Program  for  statistical  data  analysis.  (Version  2.2) 

STARPAC 

145 

Fortran  subprograms  for  statistical  data  analysis.  (Version  2.07) 

^ n is  the  number  o 

modules  classified. 
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4.1  Guidelines  for  Classifying  Modules 

When  using  the  classification  system  one  must  decide  what  objects  to  classify  and  then  how 
to  assign  classifications.  The  following  guidehnes  were  used  to  classify  software  modules  in 
the  GAMS  catalog. 

Modules  may  represent  different  types  of  objects. 

The  use  of  the  GAMS  Classification  System  is  most  straightforward  in  the  case  of  subpro- 
gram libraries,  where  one  classifies  the  individual  user-callable  subprograms.  However,  most 
statistical  software,  and  an  increasing  amount  of  mathematical  software,  comes  in  the  form 
of  stand-alone  programs  with  their  own  input  command  languages.  When  such  a program 
is  designed  to  solve  a very  restricted  set  of  problems,  then  it  is  reasonable  to  classify  it  just 
as  a subprogram  would  be.  This  is  the  case,  for  example,  with  the  programs  in  the  BMDP 
package.  In  other  cases,  a single  program  may  be  capable  of  solving  a very  wide  range  of 
problems;  interactive  statistical  analysis  systems  hke  Dataplot  and  Minitab  are  examples. 
In  such  cases  one  is  faced  with  a dilemma:  one  must  either  classify  the  program  at  a very 
high  level  (where  it  will  likely  not  be  found)  or  give  it  many  classifications  (each  of  which 
provides  httle  information).  We  have  chosen  instead  to  classify  the  major  commands  in 
the  input  language  of  these  programs.  In  this  way,  classification-based  software  advisory 
systems  have  information  available  about  the  commands  available  in  these  multi-purpose 
programs  in  the  same  way  as  user- callable  subprograms  for  the  same  problem.  This  avoids 
classification  at  too  high  a level  and  provides  catalog  readers  with  more  than  the  name  of 
the  program.  In  other  cases  we  have  chosen  to  classify  an  entire  subprogram  library  as  a 
single  unit.  This  occurs  in  the  case  of  Fortran-callable  graphics  hbraries.  We  wish  to  catalog 
the  available  graphics  libraries  in  GAMS,  but  the  classification  system  is  not  yet  sufficiently 
refined  to  meaningfully  classify  the  individual  subprograms.  In  addition,  since  users  rarely 
use  more  than  one  graphics  hbrary  at  a time,  further  refinement  is  less  critical  than  with 
mathematical  and  statistical  libraries. 

Some  modules  should  not  be  classified. 

It  is  a mistake  to  attempt  to  classify  every  software  module  in  a given  package,  even  when 
each  module  is  “user- callable”.  Some  modules,  although  called  by  the  user,  are  subsidiary 
in  nature.  Examples  of  these  are  modules  which  are  called  to  change  defaults  or  perform 
some  initialization  in  preparation  for  the  use  of  another  module.  Other  examples  are  modules 
which  evaluate  fitted  functions  or  interpolate  the  computed  solution  to  a differential  equation. 
The  detailed  documentation  of  the  module  which  solves  the  main  problem  of  interest  will 
point  to  such  subsidiary  routines.  Classifying  these  would  only  add  needless  clutter  to  a 
software  catalog. 

Modules  may  be  classified  at  more  than  one  node. 

Many  software  modules  have  multiple  purposes,  and  hence  should  be  assigned  multiple 
classes.  Also,  since  there  is  considerable  overlap  in  many  areas  of  the  classification  scheme. 
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it  may  be  difficult  to  assign  a unique  classification  in  every  case.  For  example,  nonlinear 
least  squares  approximation  (Klb)  and  nonhnear  regression  (L8b)  are  the  same  problem  seen 
from  different  points  of  view. 

It  is  interesting  to  note  that  the  assignment  of  alternate  classifications  to  modules  is 
context-sensitive,  i.e.,  it  depends  on  the  collection  of  software  being  classified.  For  example, 
consider  the  relationship  between  software  for  nonlinear  least  squares  problems  and  software 
for  nonhnear  regression.  Both  solve  the  same  basic  mathematical  problem,  but  the  latter 
software  takes  a statistical  point  of  view  (i.e.,  it  uses  the  terminology  of  statistics  and  returns 
additional  statistical  information  which  can  be  used  to  judge  goodness  of  fit).  If  the  software 
collection  is  rich  in  codes  of  both  types,  then  nonhnear  regression  codes  should  only  be 
classified  in  subtree  L8.  Users  with  a statistical  orientation  will  natural  go  down  the  L8  path 
to  find  such  codes;  assigning  L8  classes  to  nonhnear  least  squares  codes  will  only  complicate 
the  selection  process.  On  the  other  hand,  if  the  software  collection  contains  few  codes  for 
nonhnear  regression,  then  assigning  classes  in  L8  to  the  nonhnear  least  squares  codes  might 
be  the  only  way  that  naive  users  could  discover  software  appropriate  for  nonlinear  regression. 


Modules  may  be  classified  at  any  node  in  the  tree. 

We  classify  modules  at  the  lowest  level  of  the  tree  which  accurately  describes  the  problem 
solved.  This  gives  the  best  match  between  classification  and  software.  In  some  cases  the 
node  selected  in  this  manner  is  not  a leaf  of  the  tree.  This  situation  occurs  when  there  is  no 
child  node  which  adequately  describes  the  function  of  the  module.  Classification  at  non-leaf 
nodes  can  also  be  done  when  a module  solves  all  (or  at  least  most)  of  the  problems  given  at 
the  lower  levels.  This  is  not  generally  recommended  however,  since  users  will  tend  to  look 
as  far  down  into  the  tree  as  possible  in  locating  their  problems. 


4.2  Associated  Materials 

In  order  to  make  best  use  of  the  classification  system  in  the  development  of  the  GAMS 
catalog,  we  have  developed  some  associated  materials.  Among  these  are  a verbose  realization 
of  the  classification  system  and  a keyword  index  to  the  classification  system. 

The  wording  of  classes  given  in  Appendix  A is  appropriate  for  use  when  the  entire 
system  is  displayed  in  outline  form  as  it  is  there.  The  problem  descriptions  are  quite  terse 
and  are  meant  to  be  read  in  context.  In  some  cases  it  is  necessary  to  use  such  descriptions 
out  of  context,  however.  For  example,  when  modules  are  listed  in  order  of  classifications 
in  the  GAMS  catalog,  the  parent  classification  is  rarely  found  on  the  same  page  as  a given 
classification.  Similarly,  GAMSIC  displays  the  current  class  and  its  children,  prompting  the 
user  for  the  next  subclass  to  go  to.  In  each  of  these  cases  the  terse  class  description  may  not 
be  enough  for  a reader  to  determine  what  problem  the  current  class  represents.  As  a result 
of  this  we  maintain  two  separate  versions  of  the  GAMS  classification  system  which  we  refer 
to  as  terse  and  verbose.  The  terse  version  is  given  in  Appendix  A.  The  verbose  descriptions 
provide  enough  additional  wording  so  that  the  problem  represented  by  each  class  can  be 
determined  out  of  context.  The  latter  can  be  found  in  the  Modules  by  Class  section  of  [6]. 
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One  of  the  roadblocks  encountered  by  users  in  trying  to  find  software  using  the  GAMS 
Classification  System  is  that  the  terminology  used  to  describe  mathematical  and  statistical 
problems  may  not  be  familiar  to  them.  A partial  solution  developed  for  the  GAMS  catalog 
has  been  the  development  of  a keyword  index  to  the  classification  system.  Not  only  does 
this  allow  us  to  provide  pointers  into  the  classification  system  using  alternate  terminology, 
but  it  also  provides  users  a faster  method  of  getting  close  to  a desired  class  than  a linear 
search  on  paper  or  an  automated  tree  traversal.  The  index  which  we  have  developed  can  be 
found  in  [6]. 

5 Future  Development 

Since  its  creation  nearly  a decade  ago,  the  GAMS  Classification  System  has  undergone 
substantial  enhancement  and  revision.  There  is  still  much  work  to  be  done.  In  this  section 
we  describe  some  of  our  ideas  for  future  editions. 

Many  leaves  of  the  C subtree  (Elementary  and  Special  Functions)  combine  a number  of 
related  functions  in  a single  class.  For  example,  Bessel  functions  J,  Y,  Hi,  and  H2  are  all  in 
a single  class  as  are  the  Airy  functions  Ai  and  Bi.  The  amount  of  software  available  for  the 
evaluation  of  special  functions  is  steadily  increasing,  as  can  be  seen  in  the  current  GAMS 
catalog.  In  order  to  reduce  the  number  of  modules  in  these  classes  to  a more  manageable 
size  it  may  be  necessary  to  refine  many  of  the  subtrees  class  C. 

There  is  no  reasonable  place  to  classify  software  for  the  manipulation  of  piecewise  poly- 
nomials (splines)  in  the  current  scheme.  Such  software  is  currently  classified  in  subtrees 
E6  (Service  routines  for  interpolation)  or  K6  (Service  routines  for  approximation).  These 
classes  primarily  were  designed  to  contain  subsidiary  software  associated  with  interpolation 
and  approximation.  The  existence  of  these  classes  violates  our  guideline  about  not  classify- 
ing subsidiary  software  modules,  and  hence  we  intend  to  delete  these  two  subtrees.  However, 
we  must  find  a new  home  for  general-purpose  programs  for  the  manipulation  of  piecewise 
polynomials.  We  believe  that  such  a class  probably  belongs  in  C (Elementary  and  special 
functions)  in  parallel  to  class  C3  (Polynomials). 

The  amount  of  software  for  solving  problems  in  linear  algebra  is  still  increasing,  and 
many  new  problems  are  being  addressed  by  software  in  this  area.  Examples  are:  software  for 
elementary  vector  and  matrix  operations  not  listed  in  Dl  (e.g.  scalar  addition/subtraction, 
distance  between  vectors,  angles  between  vectors,  bilinear  forms),  and  software  for  the  so- 
lution of  specialized  types  of  linear  systems  (e.g.  Toeplitz,  block  tridiagonal).  Another 
problem  is  that  software  for  computing  matrix  factorizations  and  inverses  are  found  in  the 
same  classes  as  modules  for  solving  systems  of  linear  equations.  This  nearly  triples  the  num- 
ber of  modules  in  these  classes  which  are  already  among  the  most  heavily  populated  in  the 
GAMS  catalog.  Subtree  D may  need  a substantial  revision  to  alleviate  these  problems. 

The  subtree  II a (Initial- value  problems  for  ordinary  differential  equations)  violates  our 
philosophy  of  partitioning  by  problem  rather  than  by  solution  method.  The  growing  collec- 
tion of  software  for  this  problem  area  requires  us  to  find  better  ways  to  partition  it. 

Finally,  important  classes  remain  unrefined,  including  Q (Graphics)  and  0 (Symbolic 
Computation). 


The  GAMS  Problem  Classification  System 


11 


6 Conclusions 

Further  development  of  mathematical  software  advisory  systems  is  necessary  in  order  to 
ease  user  access  to  the  steadily  increasing  collection  of  reusable  mathematical  and  statistical 
software.  Tree- structured  problem-oriented  software  classification  schemes  are  one  way  for 
such  advisory  systems  to  systematically  associate  software  modules  with  the  problems  they 
solve.  Versions  of  the  Guide  to  Available  Mathematical  Software  Classification  System  have 
been  successfully  used  for  this  purpose  for  about  ten  years. 

We  seek  constructive  criticism  of  our  system,  especially  from  those  who  have  used  it  to 
classify  software.  Numerous  changes  for  the  system  are  already  being  planned;  are  seeking 
interested  parties  to  review  them.  Machine-readable  copies  of  the  system  are  available  from 
the  authors,  as  well  as  our  classifications  for  the  libraries  listed  in  Table  1. 


Disclaimer 

Certain  commercial  products  are  identified  in  this  report  in  order  to  adequately  document 
the  development  and  evaluation  of  the  GAMS  classification  system.  Identification  of  these 
products  does  not  imply  recommendation  or  endorsement  by  NIST,  nor  does  it  imply  that 
the  identified  products  are  necessarily  the  best  available  for  the  purpose. 
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Appendix  A 

GAMS  Classification  Scheme,  Version  2.0 


A.  Arithmetic,  error  analysis 

Al.  Integer 

A2.  Rational 

A3.  Real 

A3a.  Standard  precision 

A3c.  Extended  precision 

A3d.  Extended  range 

A4.  Complex 

A4a.  Standard  precision 

A4c.  Extended  precision 

A4d.  Extended  range 

A5.  Interval 

A6.  Change  of  representation 

A6a.  Type  conversion 

A6b.  Base  conversion 

A6c.  Decomposition,  construction 

A7.  Sequences  (e.g.,  convergence  acceleration) 

B.  Number  theory 

C.  Elementary  and  special  functions  (search  also  class  L5) 

Cl.  Integer-valued  functions  (e.g.,  factorial,  binomial  coefficient,  permutations,  combinations, 

floor,  ceiling) 

C2.  Powers,  roots,  reciprocals 

C3.  Polynomials 

C3a.  Orthogonal 

C3al.  Trigonometric 

C3a2.  Chebyshev,  Legendre 

C3a3.  Laguerre 

C3a4.  Hermite 

C3b.  Non-orthogonal 

C4.  Elementary  transcendental  functions 

C4a.  Trigonometric,  inverse  trigonometric 

C4b.  Exponential,  logarithmic 

C4c.  Hyperbolic,  inverse  hyperbolic 

C4d.  Integrals  of  elementary  transcendental  functions 

C5.  Exponential  and  logarithmic  integrals 

C6.  Cosine  and  sine  integrals 

C7.  Gamma 

C7a.  Gamma,  log  gamma,  reciprocal  gamma 

C7b.  Beta,  log  beta 

C7c.  Psi  function 

C7d.  Polygamma  function 

C7e.  Incomplete  gamma 

C7f.  Incomplete  beta 
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C7g. 

C8. 

C8a. 

C8b. 

C8c. 

C9. 

CIO. 

ClOa. 

ClOal. 

C10a2. 

ClOaS. 

C10a4. 

ClOaS. 

ClOb. 

ClObl. 

C10b2. 

ClObS. 

C10b4. 

ClObS. 

ClOc. 

ClOd. 

ClOe. 

Cl  Of. 

Cll. 

C12. 

C13. 

C14. 

CIS. 

C16. 

C17. 

C18. 

C19. 

D. 

Dl. 

Dla. 

Dial. 

Dla2. 

Dla3. 

Dla3a. 

Dla3b. 

Dla3c. 

Dla4. 

DlaS. 

Dla6. 

Dla7. 

Dla8. 

Dla9. 

DlalO. 

Dlall. 

Dlb. 

Dlbl. 


Riemann  zeta 
Error  functions 

Error  functions,  their  inverses,  integrals,  including  the  normal  distribution  function 
Fresnel  integrals 
Dawson’s  integral 
Legendre  functions 
Bessel  functions 
J,  Y,  Hi,  H2 

Real  argument,  integer  order 
Complex  argument,  integer  order 
Real  argument,  real  order 
Complex  argument,  real  order 
Complex  argument,  complex  order 
I,  K 

Real  argument,  integer  order 
Complex  argument,  integer  order 
Real  argument,  real  order 
Complex  argument,  real  order 
Complex  argument,  complex  order 
Kelvin  functions 
Airy  and  Scorer  functions 
Struve,  Anger,  and  Weber  functions 
Integrals  of  Bessel  functions 
Confluent  hypergeometric  functions 
Coulomb  wave  functions 
Jacobian  elliptic  functions,  theta  functions 
Elliptic  integrals 
Weierstrass  elliptic  functions 
Parabolic  cylinder  functions 
Mathieu  functions 
Spheroidal  wave  functions 
Other  special  functions 

Linear  Algebra 

Elementary  vector  and  matrix  operations 
Elementary  vector  operations 
Set  to  constant 

Minimum  and  maximum  components 
Norm 

Li  (sum  of  magnitudes) 

L2  (Euclidean  norm) 

Loo  (maximum  magnitude) 

Dot  product  (inner  product) 

Copy  or  exchange  (swap) 

Multiplication  by  scalar 

Triad  (ax  + y for  vectors  x,  y and  scalar  a) 

Elementary  rotation  (Givens  transformation) 

Elementary  reflection  (Householder  transformation) 

Convolutions 
Other  vector  operations 
Elementary  matrix  operations 

Initialize  (e.g.,  to  zero  or  identity) 
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Dlb2. 

Dlb3. 

Dlb4. 

Dlb5. 

Dlb6. 

Dlb7. 

Dlb8. 

Dlb9. 

DlblO. 

Dlbll. 

D2. 

D2a. 

D2al. 

D2a2. 

D2a2a. 

D2a3. 

D2a4. 

D2b. 

D2bl. 

D2bla. 

D2blb. 

D2b2. 

D2b2a. 

D2b4. 

D2c. 

D2cl. 

D2c2. 

D2c2a. 

D2c3. 

D2c4. 

D2d. 

D2dl. 

D2dla. 

D2dlb. 

D2d2. 

D2d2a. 

D2d4. 

D2e. 

D3. 

D3a. 

D3al. 

D3a2. 

D3a2a. 

D3a3. 

D3a4. 

D3b. 

D3bl. 

D3bla. 

D3blb. 

D3b2. 

D3b2a. 

D3b4. 


Norm 

Transpose 

Multiplication  by  vector 
Addition,  subtraction 
Multiplication 
Matrix  polynomial 
Copy 

Storage  mode  conversion 

Elementary  rotation  (Givens  transformation) 

Elementary  reflection  (Householder  transformation) 

Solution  of  systems  of  linear  equations  (including  inversion,  LU  and  related  decompositions) 

Real  nonsymmetric  matrices 
General 
Banded 

TridiagonaJ 

Triangular 

Sparse 

Real  symmetric  matrices 
General 

Indefinite 
Positive  definite 
Positive  definite  banded 
Tridiagonal 
Sparse 

Complex  non-Hermitian  matrices 
General 
Banded 

TridiagonaJ 

Triangular 

Sparse 

Complex  Hermitian  matrices 
General 
Indefinite 
Positive  definite 
Positive  definite  banded 
Tridiagonal 
Sparse 

Associated  operations  (e.g.,  matrix  reorderings) 

Determinants 

Real  nonsymmetric  matrices 
General 
Banded 

TridiagonaJ 

Triangular 

Sparse 

Real  symmetric  matrices 
General 

Indefinite 
Positive  definite 
Positive  definite  banded 
Tridiagonal 
Sparse 
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D3c. 

D3cl. 

D3c2. 

D3c2a. 

D3c3. 

D3c4. 

D3d. 

D3dl. 

D3dla. 

D3dlb. 

D3d2. 

D3d2a. 

D3d4. 

D4. 

D4a. 

D4al. 

D4a2. 

D4a3. 

D4a4. 

D4a5. 

D4a6. 

D4a7. 

D4b. 

D4bl. 

D4b2. 

D4b3. 

D4b4. 

D4b5. 

D4c. 

D4cl. 

D4cla. 

D4clb. 

D4clbl 

D4clb2 

D4clb3 

D4clc. 

D4c2. 

D4c2a. 

D4c2b. 

D4c2c. 

D4c3. 

D4c4. 

D4c5. 

D5. 

D6. 

D7. 

D7a. 

D7b. 

D7c. 

D7d. 

D8. 

D9. 


Complex  non-Hermitian  matrices 
General 
Banded 

Tridiagonal 

Triangular 

Sparse 

Complex  Hermitian  matrices 
General 
Indefinite 
Positive  definite 
Positive  definite  banded 
TridiagonaJ 
Sparse 

Eigenvalues,  eigenvectors 

Ordinary  eigenvalue  problems  [Ax  = Ax) 

Real  symmetric 

Real  nonsymmetric 

Complex  Hermitian 

Complex  non-Hermitian 

Tridiagonal 

Banded 

Sparse 

Generalized  eigenvalue  problems  (e.g.,  Ax  = XBx) 

Real  symmetric 
Real  general 
Complex  Hermitian 
Complex  general 
Banded 

Associated  operations 
Transform  problem 
Balance  matrix 
Reduce  to  compact  form 
Tridiagonal 
Hessenberg 
Other 

Standardize  problem 

Compute  eigenvalues  of  matrix  in  compact  form 
Tridiagonal 
Hessenberg 
Other 

Form  eigenvectors  from  eigenvalues 
Back  transform  eigenvectors 
Determine  Jordan  normal  form 
QR  decomposition,  Gram-Schmidt  orthogonalization 
Singular  value  decomposition 
Update  matrix  decompositions 
LU 

Cholesky 

QR 

Singular  value 

Other  matrix  equations  (e.g.,  AX  + XB  = C) 

Singular,  overdetermined  or  underdetermined  systems  of  linear  equations,  generalized  inverses 
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D9a. 

D9al. 

D9a2. 

D9a3. 

D9a4. 

D9b. 

D9bl. 

D9b2. 

D9b3. 

D9b4. 

D9c. 

E. 

El. 

Ela. 

Elb. 

Elc. 

E2. 

E2a. 

E2b. 

E3. 

E3a. 

E3al. 

E3a2. 

E3a3. 

E3b. 

E3c. 

E3d. 

F. 

FI. 

Fla. 

Flal. 

Fla2. 

Fib. 

F2. 

F3. 

G. 

Gl. 

Gla. 

Glal. 

Glala. 

Glalb. 

Glalc. 

Gla2. 

Gib. 

Glbl. 

Glbla. 

Glblb. 


Unconstrained 

Least  squares  (Z2)  solution 
Chebyshev  (Loo)  solution 
Least  absolute  value  (Li)  solution 
Other 

Constrained 

Least  squares  (L2)  solution 
Chebyshev  (Loo)  solution 
Least  absolute  value  (Li) 

Other 

Generalized  inverses 

Interpolation 

Univariate  data  (curve  fitting) 

Polynomial  splines  (piecewise  polynomials) 

Polynomials 

Other  functions  (e.g.,  rational,  trigonometric) 

Multivariate  data  (surface  fitting) 

Gridded 

Scattered 

Service  routines  for  interpolation 

Evaluation  of  fitted  functions,  including  quadrature 
Function  evaluation 
Derivative  evaluation 
Quadrature 

Grid  or  knot  generation 

Manipulation  of  basis  functions  (e.g.,  evaluation,  change  of  basis) 
Other 

Solution  of  nonlinear  equations 

Single  equation 
Polynomial 

Real  coefficients 
Complex  coefficients 
Nonpolynomial 

System  of  equations 

Service  routines  (e.g.,  check  user-supplied  derivatives) 

Optimization  (search  also  classes  K,  L8) 

Unconstrained 

Univariate 

Smooth  function 

User  provides  no  derivatives 
User  provides  first  derivatives 
User  provides  first  and  second  derivatives 
General  function  (no  smoothness  assumed) 

Multivariate 

Smooth  function 

User  provides  no  derivatives 
User  provides  first  derivatives 
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Glblc. 

Glb2. 

G2. 

G2a. 

G2al. 

G2a2. 

G2b. 

G2c. 

G2cl. 

G2c2. 

G2c3. 

G2c4. 

G2c5. 

G2c6. 

G2c7. 

G2d. 

G2dl. 

G2d2. 

G2d3. 

G2d3a. 

G2d3b. 

G2d4. 

G2e. 

G2el. 

G2e2. 

G2f. 

G2g. 

G2h. 

G2hl. 

G2hla. 

G2hlal. 

G2hla2. 

G2hla3. 

G2hlb. 

G2h2. 

G2h2a. 

G2h2al. 

G2h2a2. 

G2h2a3. 

G2h2b. 

G2h3. 

G2h3a. 

G2h3al. 

G2h3ala. 

G2h3alb. 

G2h3alc. 

G2h3a2. 

G2h3b. 

G2h3bl. 

G2h3bla. 

G2h3blb. 

G2h3blc. 


User  provides  first  and  second  derivatives 
General  function  (no  smoothness  assumed) 

Constrained 

Linear  programming 

Dense  matrix  of  constraints 
Sparse  matrix  of  constraints 
Transportation  and  assignments  problem 
Integer  programming 
Zero/one 

Covering  and  packing  problems 
Knapsack  problems 
Matching  problems 

Flouting,  scheduling,  location  problems 
Pure  integer  programming 
Mixed  integer  programming 
Network  (for  network  reliability  search  class  M) 

Shortest  path 
Minimum  spanning  tree 
Maximum  flow 

Generalized  networks 
Networks  with  side  constraints 
Test  problem  generation 
Quadratic  programming 

Positive  definite  Hessian  (i.e.,  convex  problem) 

Indefinite  Hessian 
Geometric  programming 
Dynamic  programming 
General  nonlinear  programming 
Simple  bounds 

Smooth  function 

User  provides  no  derivatives 
User  provides  first  derivatives 
User  provides  first  and  second  derivatives 
General  function  (no  smoothness  assumed) 

Linear  equality  or  inequality  constraints 
Smooth  function 

User  provides  no  derivatives 
User  provides  first  derivatives 
User  provides  first  and  second  derivatives 
General  function  (no  smoothness  assumed) 

Nonlinear  constraints 

Equality  constraints  only 

Smooth  function  and  constraints 
User  provides  no  derivatives 

User  provides  first  derivatives  of  function  and  constraints 
User  provides  first  and  second  derivatives  of  function  and  constraints 
General  function  and  constraints  (no  smoothness  assumed) 

Equality  and  inequality  constraints 
Smooth  function  and  constraints 
User  provides  no  derivatives 

User  provides  first  derivatives  of  function  and  constraints 

User  provides  first  and  second  derivatives  of  function  and  constraints 
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G2h3b2. 

G2i. 

G3. 

G4. 

G4a. 

G4b. 

G4c. 

G4d. 

G4e. 

G4f. 

H. 

HI. 

H2. 

H2a. 

H2al. 

H2ala. 

H2alal. 

H2ala2. 

H2alb. 

H2albl. 

H2alb2. 

H2a2. 

H2a2a. 

H2a2al. 

H2a2a2. 

H2a2b. 

H2a2bl. 

H2a2b2. 

H2a3. 

H2a3a. 

H2a3al. 

H2a3a2. 

H2a4. 

H2a4a. 

H2a4al. 

H2a4a2. 

H2b. 

H2bl. 

H2bla. 

H2blal. 

H2bla2. 

H2blb. 

H2blbl. 

H2blb2. 

H2b2. 

H2b2a. 

H2b2al. 

H2b2a2. 

H2b2b. 

H2b2bl. 


General  function  and  constraints  (no  smoothness  assumed) 

Global  solution  to  nonconvex  problems 
Optimal  control 
Service  routines 

Problem  input  (e.g.,  matrix  generation) 

Problem  scaling 

Check  user-supplied  derivatives 

Find  feasible  point 

Check  for  redundancy 

Other 

Differentiation,  integration 

Numerical  differentiation 

Quadrature  (numerical  evaluation  of  definite  integrals) 

One-dimensional  integrals 

Finite  interval  (general  integrand) 

Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Integrand  available  only  on  grid 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Finite  interval  (specific  or  special  type  integrand  including  weight  functions,  oscillating 

and  singular  integrands,  principal  value  integrals,  splines,  etc.) 

Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Integrand  available  only  on  grid 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Semi-infinite  interval  (including  exp  — x weight  function) 

Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Infinite  interval  (including  exp  — weight  function) 

Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 
Multidimensional  integrals 

One  or  more  hyper-rectangular  regions  (includes  iterated  integrals) 

Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Integrand  available  only  on  grid 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

n-dimensional  quadrature  on  a nonrectangular  region 
Integrand  available  via  user-defined  procedure 

Automatic  (user  need  only  specify  required  accuracy) 

Nonautomatic 

Integrand  available  only  on  grid 

Automatic  (user  need  only  specify  required  accuracy) 
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H2b2b2. 

H2c. 

I. 

11. 

Ila. 

Hal. 

Ilala. 

Ilalb. 

Hale. 

Ha2. 

Hb. 

Hbl. 

Hb2. 

Hb3. 

He. 

12. 

I2a. 

I2al. 

I2ala. 

I2alb. 

I2a2. 

I2b. 

I2bl. 

I2bla. 

I2blal. 

I2blala. 

I2blalb. 

I2bla2. 

I2bla3. 

I2ble. 

I2b2. 

I2b3. 

I2b4. 

I2b4a. 

I2b4b. 

13. 

J. 

Jl. 

Jla. 

Jlal. 

Jla2. 

Jla3. 

Jib. 

J2. 

J3. 

J4. 


Nonautomatic 

Service  routines  (e.g.,  compute  weights  and  nodes  for  quadrature  formulas) 

Differential  and  integral  equations 

Ordinary  differential  equations  (ODE’s) 

Initial  value  problems 

General,  nonstiff  or  mildly  stiff 

One-step  methods  (e.g.,  Runge-Kutta) 

Multistep  methods  (e.g.,  Adams  predictor-corrector) 

Extrapolation  methods  (e.g.,  Bulirsch-Stoer) 

Stiff  and  mixed  algebraic-  differential  equations 
Multipoint  boundary  value  problems 
Linear 
Nonlinear 

Eigenvalue  (e.g.,  Sturm-Liouville) 

Service  routines  (e.g.,  interpolation  of  solutions,  error  handling,  test  programs) 

Partial  differential  equations 

Initial  boundary  value  problems 
Parabolic 

One  spatial  dimension 
Two  or  more  spatial  dimensions 
Hyperbolic 

Elliptic  boundary  value  problems 
Linear 

Second  order 

Poisson  (Laplace)  or  Helmholtz  equation 

Rectangular  domain  (or  topologically  rectangular  in  the  coordinate  system) 
Nonrectangular  domain 
Other  separable  problems 
Nonseparable  problems 
Higher  order  equations  (e.g.,  biharmonic) 

Nonlinear 
Eigenvalue 
Service  routines 

Domain  triangulation  (search  also  class  P) 

Solution  of  discretized  elliptic  equations 
Integral  equations 

Integral  transforms 

Trigonometric  transforms  including  fast  Fourier  transforms 
One-dimensional 
Real 

Complex 

Sine  and  cosine  transforms 
Multidimensional 
Convolutions 
Laplace  transforms 
Hilbert  transforms 
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K. 

Kl. 

Kla. 

Klal. 

Klala. 

Klalal. 

Klala2. 

KlalaS. 

Klalb. 

Kla2. 

Kla2a. 

Kla2b. 

Klb. 

Klbl. 

Klbla. 

Klblal. 

Klbla2. 

KlblaS. 

Klblb. 

Klb2. 

Klb2a. 

Klb2b. 

K2. 

K3. 

K4. 

K5. 

K6. 

K6a. 

K6al. 

K6a2. 

K6a3. 

K6b. 

K6c. 

K6d. 

L. 

LI. 

Lla. 

Llal. 

Llala. 

Llalb. 

Llalc. 

Llald. 

Llale. 

Lla3. 

Lib. 

Lie. 

Llcl. 

Llclb. 

Llcld. 

Llc2. 


Approximation  (search  also  class  L8) 

Least  squares  (L2)  approximation 

Linear  least  squares  (search  also  classes  D5,  D6,  D9) 
Lfneonstrained 

Univariate  data  (curve  fitting) 

Polynomial  splines  (piecewise  polynomials) 
Polynomials 

Other  functions  (e.g.,  trigonometric,  user-specified) 
Multivariate  data  (surface  fitting) 

Constrained 

Linear  constraints 
Nonlinear  constraints 
Nonlinear  least  squares 
Unconstrained 

Smooth  functions 

User  provides  no  derivatives 
User  provides  first  derivatives 
User  provides  first  and  second  derivatives 
General  functions 
Constrained 

Linear  constraints 
Nonlinear  constraints 
Minimajc  (Loo)  approximation 
Least  absolute  value  (Li)  approximation 

Other  analytic  approximations  (e.g.,  Taylor  polynomial,  Pade) 
Smoothing 

Service  routines  for  approximation 

Evaluation  of  fitted  functions,  including  quadrature 
Function  evaluation 
Derivative  evaluation 
Quadrature 

Grid  or  knot  generation 

Manipulation  of  basis  functions  (e.g.,  evaluation,  change  of  basis) 
Other 

Statistics,  probability 

Data  summarization 

One-dimensional  data 
Raw  data 
Location 
Dispersion 
Shape 

Frequency,  cumulative  frequency 
Ties 

Grouped  data 

Two  dimensional  data  (search  also  class  Lie) 

Multi-dimensional  data 
Raw  data 

Covariance,  correlation 
Frequency,  cumulative  frequency 
Raw  data  containing  missing  values  (search  also  class  Llcl) 
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L2. 

L2a. 

L2b. 

L2c. 

L2d. 

L2e. 

L3. 

L3a. 

L3al. 

L3a2. 

L3a3. 

L3a4. 

L3a5. 

L3a6. 

L3a7. 

L3b. 

L3bl. 

L3b2. 

L3b3. 

L3b3a. 

L3b3b. 

L3b3c. 

L3b4. 

L3c. 

L3e. 

L3el. 

L3e2. 

L3e3. 

L3e3a. 

L3e3c. 

L3e3d. 

L3e4. 

L4. 

L4a. 

L4al. 

L4ala. 

L4alal. 

L4ala2. 

L4ala2b. 

L4ala2c. 

L4ala2d. 

L4ala2e. 

L4ala2f. 

L4ala2g. 

L4ala2h. 

L4ala2l. 

L4ala2n. 

L4ala2p. 

L4ala2s. 

L4ala2t. 

L4ala2u. 

L4ala2w. 


Data  manipulation 

Transform  (search  also  classes  LlOal,  N6,  and  N8) 

Tally 

Subset 

Merge  (search  also  class  N7) 

Construct  new  variables  (e.g.,  indicator  variables) 

Elementary  statistical  graphics  (search  also  class  Q) 

One-dimensional  data 
Histograms 

Frequency,  cumulative  frequency,  percentile  plots 
EDA  (e.g.,  box-plots) 

Bar  charts 
Pie  charts 

Xi  vs.  i (including  symbol  plots) 

Lag  plots  (e.g.,  plots  of  Xi  vs.  X,_i) 

Two-dimensional  data  (search  also  class  L3e) 

Histograms  (superimposed  and  bivariate) 

Frequency,  cumulative  frequency 
Scatter  diagrams 
Y vs.  X 
Symbol  plots 

Lag  plots  (i.e.,  plots  of  JV,  vs.  Y,_;  ) 

EDA 

Three-dimensional  data  (search  also  class  L3e) 

Multi-dimensional  data 
Histograms 

Frequency,  cumulative  frequency,  percentile  plots 
Scatter  diagrams 

Superimposed  Y vs.  X 
Superimposed  Xi  vs.  i 
Matrices  of  bivariate  scatter  diagrams 
EDA 

Elementary  data  analysis 
One-dimensional  data 
Raw  data 

Parametric  analysis 

Plots  of  empirical  and  theoretical  density  and  distribution  functions 
Probability  plots 
Beta,  binomial 
Cauchy,  chi-squared 
Double  exponential 
Exponential,  extreme  value 
F distribution 
Gamma,  geometric 
Halfnormal 

Lambda,  logistic,  lognormal 
Negative  binomial,  normal 
Pareto,  Poisson 
Semicircular 

t distribution,  triangular 

Uniform 

Weibull 
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L4ala3. 

L4ala3c. 

L4ala3e. 

L4ala3g. 

L4ala3l. 

L4ala3n. 

L4ala3p. 

L4ala3t. 

L4ala3w. 

L4ala4. 

L4ala4b. 

L4ala4e. 

L4ala4n. 

L4ala4p. 

L4ala4u. 

L4ala4w. 

L4ala5. 

L4ala6. 

L4ala7. 

L4alb. 

L4albl. 

L4alb2. 

L4alc. 

L4ald. 

L4a3. 

L4a4. 

L4a5. 

L4b. 

L4bl. 

L4bla. 

L4blal. 

L4bla4. 

L4blb. 

L4blc. 

L4b3. 

L4b4. 

L4b5. 

L4c. 

L4cl. 

L4cla. 

L4clb. 

L4e. 

L5. 

L5a. 

L5al. 

L5alb. 

LSalc. 

L5ald. 

L5ale. 

LSalf. 

LSalg. 

L5alh. 


Probability  plot  correlation  coefficient  plots 
Chi-squared 
Extreme  value 
Gamma,  geometric 
Lambda 
Normal 

Pareto,  Poisson 
t distribution 
Weibull 

Parameter  estimates  and  tests 
Binomial 
Extreme  value 
Normal 
Poisson 
Uniform 
Weibull 

Transformation  selection  (e.g.,  for  normality) 

Tail  and  outlier  analysis 
Tolerance  limits 
Nonparametric  analysis 

Estimates  and  tests  regarding  location  (e.g.,  median),  dispersion,  and  shape 
Density  function  estimation 
Goodness-of-fit  tests 

Analysis  of  a sequence  of  numbers  (search  also  class  LlOa) 

Grouped  and/or  censored  data 
Data  sampled  from  a finite  population 
Categorical  data 

Two  dimensional  data  (search  also  class  L4c) 

Pairwise  independent  data 
Parametric  analysis 

Plots  of  empirical  and  theoretical  density  and  distribution  functions 
Parameter  estimates  and  hypothesis  tests 
Nonparametric  analysis  (e.g.,  rank  tests) 

Goodness-of-fit  tests 
Pairwise  dependent  data 
Pairwise  dependent  grouped  data 
Data  sampled  from  a finite  population 
Multi-dimensional  data  (search  also  classes  L4b  and  L7al) 

Independent  data 

Parametric  analysis 
Nonparametric  analysis 
Multiple  multi-dimensional  data  sets 
Function  evaluation  (search  also  class  C) 

Univariate 

Cumulative  distribution  functions,  probability  density  functions 
Beta,  binomial 
Cauchy,  chi-squared 
Double  exponential 

Error  function,  exponential,  extreme  value 
F distribution 
Gamma,  general,  geometric 
Halfnormal,  hypergeometric 
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L5alk. 

L5all. 

L5aln. 

L5alp. 

L5alt. 

L5alu. 

L5alv. 

L5alw. 

L5a2. 

L5a2b. 

L5a2c. 

L5a2d. 

L5a2e. 

L5a2f. 

L5a2g. 

L5a2h. 

L5a2I. 

L5a2n. 

L5a2p. 

L5a2t. 

L5a2u. 

L5a2w. 

L5b. 

L5bl. 

L5bln. 

L5b2. 

L5b2n. 

L6. 

L6a. 

L6a2. 

L6a3. 

L6a4. 

L6a5<. 

L6a6. 

L6a7. 

L6a8. 

L6al2. 

L6al4. 

L6al6. 

L6al9. 

L6a20. 

L6a21. 

L6a22. 

L6a23» 

L6b. 

L6b3. 

L6b5. 

L6bl2. 

L6bl3. 

L6bl4. 

L6bl5. 

L6b21. 


Kendall  F statistic,  Kolmogorov-Smirnov 

Lambda,  logistic,  lognormal 

Negative  binomial,  normal 

Pareto,  Poisson 

t distribution 

U niform 

Von  Mises 

Weibull 

Inverse  distribution  functions,  sparsity  functions 
Beta,  binomial 
Cauchy,  chi-squared 
Double  exponential 

Error  function,  exponential,  extreme  value 
F distribution 
Gamma,  general,  geometric 
Halfnormal 

Lambda,  logistic,  lognormal 

Negative  binomial,  normal,  normal  order  statistics 
Pareto,  Poisson 
t distribution 
Uniform 
Weibull 
Multivariate 

Cumulative  multivariate  distribution  functions,  probability  density  functions 
Normal 

Inverse  cumulative  distribution  functions 
Normal 

Random  number  generation 
Univariate 

Beta,  binomial,  Boolean 
Cauchy,  chi-squared 
Double  exponential 
Exponential,  extreme  value 
F distribution 

Gamma,  general  (continuous,  discrete),  geometric 
Halfnormal,  hypergeometric 
Lambda,  logistic,  lognormal 

Negative  binomial,  normal,  normal  order  statistics 
Pareto,  Pascal,  permutations,  Poisson 
Samples,  stable  distribution 
t distribution,  time  series,  triangular 
Uniform  (continuous,  discrete),  uniform  order  statistics 
Von  Mises 
Weibull 
Multivariate 

Contingency  table,  correlation  matrix 
Experimental  designs 

Linear  Li  (least  absolute  value)  approximation 

Multinomial 

Normal 

Orthogonal  matrix 
Uniform 
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L6c. 

L7. 

L7a. 

L7al. 

L7a2. 

L7b. 

L7c. 

L7d. 

L7dl. 

L7d2. 

L7d3. 

L7e. 

L7f. 

L7g. 

L8. 

L8a. 

L8al. 

L8ala. 

L8alal. 

L8ala2. 

L8ald. 

L8a2. 

L8a3. 

L8a4. 

L8b. 

L8bl. 

L8bla. 

L8blb. 

L8blbl 

L8blb2 

L8blc. 

L8bld. 

L8c. 

L8cl. 

L8cla. 

L8clal. 

L8cla2. 

L8cla3. 

L8clb. 

L8clbl. 

L8clb2. 

L8clc. 

L8cld. 

L8c2. 

L8c3. 

L8c4. 

L8c5. 

L8c6. 

L8d. 

L8e. 

L8el. 

L8ela. 


Service  routines  (e.g.,  seed) 

Analysis  of  variance  (including  analysis  of  covariance) 

One-way 

Parametric 

Nonparametric 

Two-way  (search  also  class  Lid) 

Three-way  (e.g.,  Latin  squares)  (search  also  class  Lid) 
Multi-way 

Balanced  complete  data  (e.g.,  factorial  designs) 

Balanced  incomplete  data 

General  linear  models  (unbalanced  data) 

Multivariate 

Generate  experimental  designs 
Service  routines 

Regression  (search  also  classes  D5,  D6,  D9,  G,  K) 

Simple  linear  (i.e.,  y = feo  + h\x)  (search  also  class  L8h) 

Ordinary  least  squares 
Parameter  estimation 
Unweighted  data 
Weighted  data 

Inference  (e.g.,  calibration)  (search  also  class  LSala) 

Lp  for  p different  from  2 (e.g.,  least  absolute  value,  minimax) 
Robust 

Errors  in  variables 

Polynomial  (e.g.,  y = 6o  -I-  bix  -I-  b2X^)  (search  also  class  L8c) 
Ordinary  least  squares 
Degree  determination 
Parameter  estimation 

Not  using  orthogonal  polynomials 
Using  orthogonal  polynomials 
Analysis  (search  also  class  L8hlb) 

Inference  (search  also  class  L8blb) 

Multiple  linear  (i.e.,  y = 6o  + biX\  -f  . . . -f  bpXp) 

Ordinary  least  squares 
Variable  selection 
Using  raw  data 

Using  correlation  or  covariance  data 
Using  other  data 

Parameter  estimation  (search  also  class  L8cla) 

Using  raw  data 
Using  correlation  data 

Analysis  (search  also  classes  L8cla  and  L8clb) 

Inference  (search  also  classes  L8cla  and  L8clb) 

Several  regressions 
Lp  for  p different  from  2 
Robust 

Measurement  error  models 
Models  based  on  ranks 
Polynomial  in  several  variables 
Nonlinear  (i.e.,  y = F(X,b))  (search  also  class  L8h) 

Ordinary  least  squares 
Variable  selection 
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L8elb. 

LSelbl. 

L8elb2. 

L8elb3. 

L8elb4. 

L802* 

L8e5. 

L8f. 

L8g. 

L8h. 

L8i. 

L9. 

L9a. 

L9b. 

L9c. 

L9d. 

LIO. 

LlOa. 

LlOal. 

LlOala. 

LlOalb. 

LlOalc. 

LlOalcl. 

Ll0alc2. 

LlOalcS. 

L10alc4. 

LlOald. 

L10a2. 

L10a2a. 

L10a2al. 

L10a2a2. 

L10a2b. 

L10a2c. 

L10a2cl. 

L10a2c2. 

L10a2d. 

L10a2dl. 

L10a2d2. 

L10a2d3. 

Ll0a2e. 

L10a2f. 

L10a3. 

L10a3a. 

Ll0a3al<. 

L10a3a2<. 

L10a3a3. 

L10a3a4. 

Ll0a3a5. 

L10a3a6. 

L10a3b. 

LlOb. 

Ll0b2. 


Parameter  estimation  (search  also  class  LSela) 

Unweighted  data,  user  provides  no  derivatives 
Unweighted  data,  user  provides  derivatives 
Weighted  data,  user  provides  no  derivatives 
Weighted  data,  user  provides  derivatives 

Ridge 

Measurement  error  models 
Simultaneous  (i.e.,  Y = Xb) 

Spline  (i.e.,  piecewise  polynomial) 

EDA  (e.g.,  smoothing) 

Service  routines  (e.g.,  matrix  manipulation  for  variable  selection) 

Categorical  data  anaJysis 
2-by-2  tables 

Two-way  tables  (search  also  class  L9d) 

Log-linear  model 

EDA  (e.g.,  median  polish) 

Time  series  analysis  (search  also  class  J) 

Univariate  (search  also  classes  L3a6  and  LSal) 

Transformations 

Elementary  (search  also  class  L2a) 

Stationarity  (search  also  class  LSal) 

Filters  (search  also  class  K5) 

Difference 

Symmetric  linear  (e.g.,  moving  averages) 

Autoregressive  linear 
Other 
Taper 

Time  domain  analysis 
Summary  statistics 

Autocorrelations  and  autocovariances 
Partial  autocorrelations 

Stationarity  analysis  (search  also  class  L10a2a) 

Autoregressive  models 
Model  identification 
Parameter  estimation 

ARMA  and  ARIMA  models  (including  Box-Jenkins  methods) 

Model  identification 
Parameter  estimation 
Forecasting 

State-space  analysis  (e.g.,  Kalman  filtering) 

Analysis  of  a locally  stationary  series 
Frequency  domain  analysis  (search  also  class  Jl) 

Spectral  analysis 
Pilot  analysis 
Periodogram  analysis 

Spectrum  estimation  using  the  periodogram 

Spectrum  estimation  using  the  Fourier  transform  of  the  autocorrelation  function 
Spectrum  estimation  using  autoregressive  models 
Spectral  windows 
Complex  demodulation 

Two  time  series  (search  also  classes  L3b3c,  LlOc,  and  LlOd) 

Time  domain  analysis 
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L10b2a. 

L10b2b. 

L10b3. 

LlObSa. 

L10b3a2. 

L10b3a3. 

L10b3a4. 

L10b3a6. 

LlOc. 

LlOd. 

Lll. 

L12. 

L13. 

L13a. 

L13b. 

L13c. 

L14. 

L14a. 

L14al. 

L14ala. 

L14alal. 

L14ala2. 

L14ala3. 

L14ala4. 

L14alb. 

L14a2. 

L14b. 

L14c. 

L14d. 

L15. 

L16. 

L17. 

M. 

Ml. 

Mia. 

Mlb. 

M2. 

M3. 

M3a. 

M3b. 

M4. 

N. 

Nl. 

N2. 

N3. 

N4. 

N5. 

N5a. 


Summary  statistics  (e.g.,  cross-correlations) 

Transfer  function  models 
Frequency  domain  analysis  (search  also  class  Jl) 

Cross-spectral  analysis 

Cross-periodogram  analysis 

Cross-spectrum  estimation  using  the  cross-periodogram 

Cross-spectrum  estimation  using  the  Fourier  transform  of  the  cross-correlation 
or  cross-covariance  function 
Spectral  functions 

Multivariate  time  series  (search  also  classes  Jl,  L3e3  and  LI  Ob) 

Two  multi-channel  time  series 
Correlation  analysis  (search  also  classes  L4  and  L13c) 

Discriminant  analysis 
Covariance  structure  models 
Factor  analysis 

Principal  components  analysis 
Canonical  correlation 
Cluster  analysis 
One-way 

Unconstrained 

Nested 

Joining  (e.g.,  single  link) 

Divisive 

Switching 

Predict  missing  values 
Non-nested  (e.g.,  K means) 

Constrained 

Two-way 

Display 

Service  routines  (e.g.,  compute  distance  matrix) 

Life  testing,  survival  analysis 
Multidimensional  scaling 
Statistical  data  sets 

Simulation,  stochastic  modeling  (search  also  classes  L6  and  LlO) 

Simulation 

Discrete 

Continuous  (Markov  models) 

Queueing 

Reliability 

Quality  control 
Electrical  network 
Project  optimization  (e.g.,  PERT) 

Data  handling  (search  also  class  L2) 

Input,  output 
Bit  manipulation 
Character  manipulation 

Storage  management  (e.g.,  stacks,  heaps,  trees) 

Searching 

Extreme  value 
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N5b. 

N5c. 

N6. 

N6a. 

N6al. 

N6ala. 

N6alb. 

N6alc. 

N6a2. 

N6a2a. 

N6a2b. 

N6a2c. 

N6b. 

N7. 

N8. 

O. 


P. 


Q. 


R. 

Rl. 

R2. 

R3. 

R3a. 

R3b 

R3c. 

R4. 

S. 

51. 

52. 

53. 

Z. 


Boisvert,  Howe  and  Kahaner 


Insertion  position 
On  a key 
Sorting 
Internal 

Passive  (i.e.  construct  pointer  array,  rank) 

Integer 

Real 

Character 

Active 

Integer 

Real 

Character 

External 

Merging 

Pernauting 

Symbolic  computation 

Computational  geometry  (search  also  classes  G and  Q) 


Graphics  (search  also  class  L3) 

Service  routines 

Machine-dependent  constants 

Error  checking  (e.g.,  check  monotonicity) 

Error  handling 

Set  criteria  for  fatal  errors 
Set  unit  number  for  error  messages 
Other  utilities 
Documentation  retrieval 

Software  development  tools 

Program  transformation  tools 
Static  program  analysis  tools 
Dynamic  program  analysis  tools 

Other 
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Appendix  B 

Changes  from  Version  1.2 

Here  we  summarize  how  the  GAMS  Classification  System  has  changed  from  version  1.2  which  was 
published  in  [3].  A number  of  additional  cosmetic  changes  were  made  to  the  text  of  the  system;  these  are 
not  enumerated  here. 

• Subtree  AS  (Real  arithmetic) 

A3a  changed  from  Real  to  Standard  precision 

• Subtree  (Complex  arithmetic) 

A4a  changed  from  Real  to  Standard  precision 

• Dlall  (Other  vector  operations) 

New  subclass. 

• D9  (Singular,  overdetermined  or  underdetermined  systems  of  linear  equations,  generalized  inverses) 
Tree  refined,  11  new  subclasses  added. 

• E3  (Service  routines  for  interpolation) 

Tree  refined,  14  new  subclasses  added. 

• Subtree  FI  (Single  nonlinear  equations) 

Revised  to  remove  distinction  between  smooth  and  nonsmooth  functions. 

• Subtree  F2  (System  of  nonlinear  equations) 

All  subclasses  deleted,  removing  the  distinction  between  smooth  and  nonsmooth  functions. 

• K6  (Service  routines  for  approximation) 

Tree  refined,  14  new  subclasses  added. 

• Subtree  L (Statistics  and  Probability) 

Substantially  revised.  Classes  LI,  L3,  L4,  L7,  and  L14  were  revised  to  standardize  the  first  level  of 
subclasses  (e.g.,  Lla,  Lib,  . . .)  as  the  dimension  of  the  data;  class  L8  was  revised  so  that  the  first  level 
of  its  subclasses  were  functional  form.  Classes  LI,  L2,  L3,  L4,  L7,  L8,  and  LIO  were  revised  to  better 
reflect  available  software.  Probability  plots  were  moved  from  L3  to  L4.  Class  L13  was  completely 
revised.  Classes  L16  and  L17  were  added. 

• Subtree  N6alb  (Internal  sorting  of  real  data) 

All  subclasses  deleted,  removing  the  distinction  between  single  and  double  precision  data. 

• Subtree  N6a2b  (Internal  sorting  of  real  data) 

All  subclasses  deleted,  removing  the  distinction  between  single  and  double  precision  data. 

• Subtree  P (Computational  Geometry) 

All  subclasses  deleted. 

• Subtree  Q (Graphics) 

Subclass  Ql  (Line  printer  graphics)  deleted.  Q has  no  subclasses  in  the  revised  system. 

In  addition,  the  text  of  the  following  classes  was  revised  in  order  to  clarify  their  purposes:  Cl,  H2bl,  H2b2, 
Jl,  Jla3. 

Classes  A3,  A4,  N6alb  and  N6a2b  all  were  changed  to  remove  classes  which  referred  to  double  precision 
to  insure  that  both  single  and  double  precision  versions  of  programs  could  always  have  the  same  classification. 

Classes  D9,  E3,  and  K6  were  refined  because  of  a wealth  of  software  now  available  for  these  problems. 
The  subtrees  E3  and  K6  are  exactly  parallel;  they  provide  homes  for  low-level  routines  for  manipulating 
spline  basis  functions. 

The  subtrees  P and  Q were  both  trimmed  so  that  they  would  be  unrefined.  In  the  case  of  P we  felt  that 
not  enough  software  was  currently  available  to  adequately  define  the  subject  area.  In  the  case  of  Q we  did 
not  feel  that  we  possessed  the  expertise  to  adequately  refine  it. 


. A3b  (Double  precision)  removed. 
. A4b  (Double  precision)  removed. 
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